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PsJ ■ Abstract 



We study the long-time behavior of global strong solutions to a hydrodynamic 
^t' ^ I system for nonhomogeneous incompressible nematic liquid crystal flows driven by two 

types of external forces in a smooth bounded domain in . For arbitrary large regular 
, initial data with the initial density being away from vacuum, we prove the decay 

of the velocity field for both cases. Furthermore, for the case with asymptotically 
p ^ '■ autonomous external force, we can prove the convergence of the density function and 

, the director vector as time goes to infinity. Estimates on convergence rate are also 

provided. 
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> ! 1 Introduction 
(N 
^ ■ 

1/-^ ■ Liquid crystals are substances that exhibit a phase of matter that has properties between 

. those of a conventional liquid, and those of a solid crystal [S]. The hydrodynamic theory of 

^ i liquid crystals due to Ericken and Leslie was developed around 1960's [51ll5pi6| . Since then, 

I the mathematical theory is still progressing and the study of the full Ericksen-Leslie model 

presents relevant mathematical difficulties. We consider the following hydrodynamical 
model for the flow of nematic liquid crystals (cf. [T7j) 

h: pt + v-vp = 0, (1.1) 



p{vt + v-Vv)-i'Av + VP = -XV ■ (VdQVd) + pg, (1.2) 
V-v = 0, (1.3) 
dt + vVd = 7(Ad-/(d)), (1.4) 

in X where fl C (n = 2,3) is assumed to be a bounded domain with smooth 
boundary F. System ()l.ip - ()1.4p is subject to the Dirichlet boundary conditions: 

v{x,t) = 0, d{x,t) = do{x), for (j;,t) G r X M+. (1.5) 
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and the initial conditions 

p\t=o = Poix), ^1^=0 = ^0(2;) with V- 170 = 0, d\t=o = do{x), for x G 0. (1.6) 

In the above system, p is the density of the material, v is the velocity field of the flow and 
d represents the averaged macroscopic/continuum molecular orientation in R". P{x,t) 
is a scalar function representing the pressure (including both the hydrostatic and the 
induced elastic part from the orientation field), g stands for the external body force. The 
positive constants A and 7 stand for viscosity, the competition between kinetic energy 
and potential energy, and macroscopic elastic relaxation time (Debroah number) for the 
molecular orientation field. We assume that f{d) = VF{d) for some smooth bounded 
function F : M" —?■ M. Vd Vd denotes the n x n matrix whose (i, j)-th entry is given by 
Vjd • Vjd, for 1 < i, j < n. 

System (|l.ip - (jl.4p is a nonhomogeneous version of the following simplified system 
introduced in [17, 20j that models the incompressible fiow of nematic liquid crystals with 
varying director lengths: 

vt + v-Vv-i^Av + VP = -AV-(Vd0Vd), (1.7) 
V-v = 0, (1.8) 
dt + vVd = j{Ad-f{d)). (1.9) 

System (jl.7p - p.9p keeps the important mathematical structure as well as some of the 
essential features of the original Ericksen-Leslie system. The Ginzburg-Landau approxi- 
mation 

f(d) = ir(\d\^ - l)d, with its antiderivative F(d) = -^(\d\^ - l)^ (1.10) 

was introduced in order to relax the nonlinear constraint \d\ = 1. The system ()1.7p - 
(jl.Op has been studied in a series of work not only theoretically [20l[2Tl[2H[33] but also 
numerically |23y25] (see also [11] for the case f{d) = 0). For mathematical results on the 
system under constraint \d\ = 1, we refer to [T9l|22l[32] and the references cited therein. 
In [20], the authors proved the existence of global weak solutions to system ()1.7p - ()1.9p 
with Dirichlet boundary conditions by a semi-Galerkin method. Global existence and 
uniqueness of classical solutions to the same system was proved for n = 2 or n = 3 under 
large viscosity assumption (see also [H] for the case with no-slip boundary conditions). 
Long-time behavior of global classical solutions to system ()1.7p - ()1.9p was studied in [201 
[27l[33]. In particular, convergence of global classical solutions to single steady states as 
time goes to infinity was obtained in [27l|33j. As far as the density-dependent case is 
concerned, the authors in [101126) proved existence of global weak solutions of the problem 
(jl.ip ~ (jl.6p without assuming the positive lower bound for the initial density. The basic 
idea of their proof is to introduce a viscous term eAp in the transport equation (jl.ip . and 
then pass to the limit as e — ?> 0. In the recent work [3J, instead of introducing the viscosity 
term in (jl.ip . the authors provided an alternative proof for the existence of global weak 
solutions to system (jl.ip - (|1.6p under a stronger assumption (jl.lip (i.e., the initial density 
is positive and bounded). Regularity properties of weak solutions to system ()l.ip - ()1.6p 
was proved by using Ladyzhenskaya type energy estimates for the approximate solutions 
constructed within a proper Galerkin scheme, provided that the initial data are regular 
and satisfy assumptions ()l.lip - ()1.12p . 
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We note that the external force pg is supposed to be vanishing in the above mentioned 
work. In this paper, we focus on the two-dimensional case n = 2 and extend the results on 
long-time behavior of global classical solutions in [201l33j to the nonhomogeneous system 
(jl.ip ~ (jl.4p with non- vanishing external forces. Two types of external forces will be treated 
in the following text: 

(Fl) g is a time-independent potential field, namely, 

g = Vi;^, for some scalar function (t>{x) S H^(Q). 

(F2) g depends on time and satisfies the following integrability conditions: 

g G L\o, +00; gt G L\o, +00; L\n)). 

The main results of this paper are as follows: 

Theorem 1.1. Suppose that r € (2, +00), the external force g satisfies either (Fl) or 
(F2). For any initial data po G W^'^'iQ), vq G 11^(0) n V and do G H^(il) that satisfy 

0<P<Po{x)<p, VxGfi, (1.11) 
|do(2;)|<l, VxgJ^, (1.12) 

where p and p are positive constants, problem ()l.ip - (|1.6p admits a unique global strong 
solution {p, V, d) such that for any T > 

pGC([o,r],w^i'^(o)), 

V G C([0,r];H2(S7) nF) nL2(o,T;H=^), vt G L'^{^,T;V), 

d G C([0, T]; H3(17)) n L2(0, T; H^), dt G C([0, T]-^l{n)) n L2(0, T; H^), 

< p< /o(x,t) < p, [d(x,t)|<l, V(x,t) G X [0,r]. 

Theorem 1.2. Suppose that the assumptions of Theorem \l.l\ are satisfied with g fulfilling 
(Fl). The global strong solution to problem (jl.ip - (|1.6p has the following property 

hm {\\v{t)\\^.+\\vt\\ + \\dt\\^.)=Q. (1.13) 

t— J- + 00 

For any unbounded sequence {ti}, there is a subsequence {t'-} /• +00 such that 

\\p{t'i) - PocWli ^ 0, as t- ^ +00, g G (1, +00), (1.14) 
\\d{t'i) - dooWm ^ 0, as t- ^ +00, (1-15) 

where poo is a certain function that belongs to L'^ and d^o is a solution to the following 
nonlinear elliptic boundary value problem: 

I -Adoo + /(doo) = 0, xefl, {116) 
\ doo = do{x), X eV. 

Theorem 1.3. Suppose that the assumptions of Theorem \l.l\ are satisfied with g fulfilling 
(F2). If in addition, g satisfies 

r+oo 

sup(l + t)^+^ / ||g(r)f dr < +00, for some £, > 0, (1.17) 
t>o Jt 
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then the global strong solution to problem. ()l.ip - ()1.6p has the following property 

lim {\\p{t) - p^Wli + \\vit)\\ui + \\vt\\ + WdtWui + Mt) - doollna) = 0, (1.18) 

where poo is certain function in L'^ [q € (1, +oo)) and d^o is a solution of ()1.16p . Moreover, 
there exists a positive constant C depending on vo,do,iy,r],^l, p, p,doo, such that 

Mt)-Poo\\H-^ + Mt)\\m + \\d{t)-doo\\n^ <C{l + t)-^, yt>0, (1.19) 
where k = min | ^^^^ , || with 9 G (O, ^) being a constant depending on /, doo- 

The question on the uniqueness of asymptotic hmit for liquid crystal system p.7p - p.9p 
was raised in |20j. A positive answer was given in [33j such that for any global classical 
solution {v, d) of the homogeneous system, the velocity field will decay to zero and the 
director vector will converge to a steady state that is a solution to the stationary problem 
(jl.lGp (cf. |27) for some generalization). Decay of the velocity field can be obtained by 
exploring the dissipative nature of the problem and by proper energy estimates, while 
convergence of the director vector is a nontrivial problem, because the structure of the 
equilibria set (namely, the set of solutions to (|1.16p ) can be quite complicated. The proof 
in |33j relies on the so-called Lojasiewicz-Simon approach [29], which turns out to be a 
useful method to study convergence of global solutions to equilibria for nonlinear evolution 
equations (see e.g., [21[S1IH1[I21[IS] and the references cited therein). One advantage of the 
approach is that, one can obtain the convergence result without studying the structure of 
the set of equilibria, which is usually difficult when the spacial dimension is larger than 
one. The nonhomogeneous problem ()l.ip - ()1.6p under consideration is much more involved 
than the homogeneous case. Main difficulties come from those nonlinear couplings between 
the three equations for density, velocity and director in terms of convection, extra stress 
term as well as the external force. Under both assumptions (Fl) and (F2), we are able to 
derive certain (dissipative) basic energy inequalities for problem ()l.ip - ()1.6p and further- 
more, some specific higher-order differential inequalities in the sprit of [20j, which not only 
provide uniform-in-time estimates for the global strong solutions but also yield decay prop- 
erty of the velocity. Our results on the decay of velocity fields under both types of external 
forces imply that the dissipations from the viscosity and the relaxation effect in (|1.4p are 
strong enough to compensate the effects of external forces and the density fluctuation as 
well as the interactions between the fluid and the liquid crystal molecules such that the 
flow will slow down as time goes to inflnity. This extends the result on density-dependent 
incompressible Navier-Stokes equations driven by a time-independent external force (cf. 
(Fl)) on bounded domains in 2D |34j . For the asymptotically autonomous external force 
(cf. (F2)), we are able to apply the Lojasiewicz-Simon approach to prove convergence of 
the director vector and thus generalize the previous results in [20y33j for the homogeneous 
liquid crystal system. In this case, we can also obtain L^-integrability of the velocity 
fleld, which together with the transport equation (jl.ip yields convergence of the density 
function. Besides, by the Lojasiewicz-Simon approach we can derive some explicit decay 
rates of the density, velocity fleld and director vector. We remark that in the case of 
time-independent external force, the Lojasiewicz-Simon method seems fail to apply, thus 
we are only able to show certain sequential convergence of the density function and the 
director vector. Our results still hold in three-dimensional case provided that bounded 
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global strong solutions of problem ()l.ip - ()1.6p can be obtained (this could be verified, for 
instance, if the initial data and the external forces are sufficiently small). 

The remaining part of this paper is organized as follows. In Section 2, we introduce the 
functional settings, some preliminary results as well as some technical lemmas. In Section 
3, we derive some dissipative energy inequalities and some specific higher-order differential 
inequalities for both types of external forces, which enable us to obtain uniform a priori 
estimates and conclude the existence and uniqueness of global strong solutions to problem 
p.l|) - (|1.6|) (cf. Theorem II. ip . Section 4 is devoted to the proof of our main results on the 
long-time behavior of global strong solutions (cf. Theorems [L2l [L3]) . 

2 Preliminaries 

As usual, LP{Q) and W'''^{Q) stand for the Lebesgue and the Sobolev spaces of real 
valued functions, with the convention that H^{Q) = W^''^{^). The spaces of vector- valued 
functions are denoted by bold letters, correspondingly. Without any further specification, 
II • II stands for the norm in L^($7) or L^(f]). We shall denote by C the genetic constants 
depending on A, 7, v, rj, $7 and the initial data. Special dependence will be pointed out 
explicitly in the text if necessary. 

For 1 < q < +00, X'^ denotes the space that is a completion in L'' of the set of 
solenoidal vector fields with coefficients that belong to V = Cg°(r2) n : V • u = 0}: 

X'i = {v L''(J7) : V • u = in J], 7; • n = on T}, 

where n is the unit outer normal to the boundary (cf. ^). In particular, we denote 

H = X'^ = the closure of V in V'iVt), V = the closure of V in Hj(n). 

It is well-known that any vector-field with coefficients in has a Helmholtz decomposition 
(cf. [28]). Denote by Pg : JJ^{Q) — > X'^ the projector from L"^ to X^, which is a bounded 
operator. We can define the Stokes operator Aq = Pq(— A) with domain D{Aq) = W^'^ n 
Wo'« n XI. The following result holds (cf. e.g., g]): 

Lemma 2.1. Let Q be a bounded domain ofM? with smooth boundary and d{i^) be the 
diameter ofQ. Then following results hold true: 

(i) For any v G W^'"? PI Wq''', there is a constant C = C{q, d{^)) such that 

||v||w2.9 := ||V^u||l'! +d{Qy'^\\Vv\\-L<i + d{Q)~'^\\v\\-L, < C||V\||l9 

(a) For 1 < q < +00, / € L'', the Stokes problem 

-Av + VP = f, in n, v\r = 0, 

has a unique solution {v,P) in D{Aq) x W^'"^. There exists a constant C = C{q,d{^)) 
such that 

||V\||l« + ||VP||l. <C||/||l.. 

Concerning the transport equation for the density function, we have the following 
result that can be found e.g., in [l] 
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Lemma 2.2. Let v € L^(0,T;Lip) be a solenoidal vector field such that v ■ n = on F. 
For any po G W^''^ with q G [l,+cx3], the equation 

pt + V -Vp = 0, p\t=o = po{x), 

admits a unique solution p £ L~(0, T; n C{[0,T];nr<ooW^''') if q 

p G C{[0,T];W^''^) if I < q < +oo. Besides, the following estimate holds 

IIpWIIh/m <e^«"^''^"^"^°°'llPo||H/M, VtG [o,r]. 

If Po G LP for some p G [l,+oo], then \\p{t)\\LP = \\po\\lp for t G [0,T]. 

An essential characteristic of the director equation (for any given given velocity) is the 
following weak maximum principle (cf. [20] ) : 

Lemma 2.3. Suppose v G L°°{0,T;H) n L'^{0,T;V), do G U.^{n) with do G Hi(r) and 

\do\ < 1- If d £ L°°(0, T; H"*^) fl -L^(0,T;H^) is the weak solution of the initial boundary 
value problem 

dt + V ■ Vd = ^{Ad — f{d)), a.e. inCl, 
d\r = do{x), gTx (0,T) 

d\t=o = doix), 

then \d{x,t)\ < 1, a.e. in Q x (0, T). 

Finally, we report some inequalities that will be frequently used in the subsequent 
proof. 

Lemma 2.4. Let Q be a bounded domain o/M^ with smooth boundary. 

(i) Ladyzhenskaya inequality. ||f||L4 < C'|| Vf || G HQ(il). 

(a) Agmon inequality. ||/|||oo < C'||/||//2 V/ G II'^{Q). 

(Hi) Poincare inequality. For 1 < q < +oo, \\v\\i,q < C\\Vv\\-L,q, G WQ''^(il). 

(iv) Sobolev embeddings. For 1 < q < +oo, the embedding ^ L"? is compact. Besides, 

II/IIl- <c||/|[h/i,„ yfew^^i{n). 

The constant C in the above inequalities may depend on Vt and q. 

3 Global strong solutions 

In order to prove the existence of strong solutions, we can first construct a sequence of 
approximate solutions {pm,Vm,dm) within a semi-Galerkin scheme as in [3l Section 5] (we 
also refer to [20j for the case of homogeneous liquid crystal flow and [IJ for the density- 
dependent Navier-Stokes equation with external forces). To prove the convergence of the 
approximate solutions, we only need to derive some a priori estimates for them. Due 
to the positivity condition (|l.lip . the Galerkin approximation in does not rely on the 
introduction of a viscosity term in the transport equation like in [TOIISB] and thus it 
can be used to establish (higher-order) Ladyzhenskaya type energy estimates. Since the 
calculations for the approximate solutions are (formally) identical to that as we work with 
smooth solutions, in what follows, we simply perform calculations for smooth solutions to 
problem I^J^-I^M). 



+00 and 



(2.20) 
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3.1 Dissipative basic energy inequalities and lower-order energy esti- 
mates 

The total energy of problem ()l.ip - (jl.6p is defined as follows: 

£{t) = - [ p{t)\v{t)\^dx + -\\Vd{t)f + X [ F{d{t))dx. (3.1) 



In analogy to the constant density case (cf. [20j) or the nonhomogeneous system without 
external force (cf. [inil26]), our system ()l.ip - (jl.6p still has the following basic energy 
inequalities, which reflect the energy dissipation of the liquid crystal flow . 

Lemma 3.1 (Basic energy inequalities). Let {p,v,d) be a smooth solution of problem 
([TTIl - ffLG]) onnx [0,T] = Qt (0 < T < +oo). 
(i) If g satisfies (Fl), it holds 

— ^(t) + z^||Vt;(t)f + A7||Ad(i) - /(d(t))||2 = 0, < t < T, (3.2) 
dt 

where 

S{t) = 8{t) - / p(l)dx. 



(a) If g satisfies (F2), it holds 

j^£{t) + '^\\Vv{t)f + A7||Ad(t) - f{d{t))f < ^||gf , 0<t<T. (3.3) 

Proof Multiplying (fTH), ([L2]) and ([131) by l\v\^, v and X{-Ad + f{d)), respectively, 
integrating over we infer from the boundary conditions (jl.6p and integration by parts 
that ^ ^ ^ 

- / pt\v\'^dx = — / \/ ■ {pv)\vf dx = - / pv ■ 'V\v\'^dx, (3-4) 

If d , .0 , 1,2 
— / P-r\'v\ dx + vllS/vW 

2.L' dr ' " " 



~2 



- / pv-V\v\^dx+ [ pg ■ vdx — \ f {Ad ■ \7d) ■ vdx, (3.5) 
2 Jn J 



^ QllVdf + F{d)dx^ + A7II - A/ + /((i)f + xj{vVd)- Addx = 

where we have used the facts (cf. [20]) 

1 



(3.6) 



V- (VdoVd) = ^V|Vdp + Ad- Vd, (3.7) 

/ VP- vdx = / V\Vd\^-vdx= / v ■VF{d)dx = 0. (3.8) 
Jn Jn Jn 

Adding ()3.4p - ()3.6p together, we can see that 



j^£{t) + i^\\Vv{t)f + Xj\\Ad{t) - fid{t))f = J^PS- vdx. (3.9) 

If g satisfies (Fl), using the idea in [31], we multiply the transport equation (jl.ip by — </> 
and integrate over $7 to get 

ptcpdx = / V ■ {pv)(j)dx = — 1 pV(j)-vdx. (3.10) 
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Adding ()3.9p with ()3.10p and noticing that (j) is independent of time, we arrive at our 
conclusion ()3.2p . On the other hand, If g satisfies (F2), using the Holder inequality and 
Poincare inequality, we infer that 



pg • vdx 



n 



< llp||L-llg|lll^ll < Cpp\\g\\\Nv\\ < ol|Vt;f + ^^||gf , (3.11) 



2" " 2iy 

which together with (|3.9p yields ()3.3p . □ 
Proposition 3.1. Under the assumptions of Theorem \l.l{ the following estimates hold 

0< P< p{x, t) <p, Vt > 0, (3.12) 

\d{x,t)\<l, Vt>0, (3.13) 

\\v{t)\\ + \\d{t)\\-tii < C, Vt>0, (3.14) 

iy\\Vv{t)f + \-f\\Ad{t) - f{d{t))f)dt < C, (3.15) 



where C is a constant depending on \\vo\\, \\do\\iii, r], p, p, VL and also \\4>\\m (under (Fl)) 
or ||g||L2(o,+oo;L2) (under (F2)). 

Proof. ()3.12p easily follows from the characteristics method while (|3.13p is a conse- 
quence of the weak maximum principle for the director equation (see Lemma 12. 3p . If g 
satisfies (Fl), since 

< < +00, (3.16) 





■pdx 


Jn 





we can see that 8{t) > —\Q,\p\\(l)\\i^i for all t > 0. The required uniform estimates follow 
from (|3.2p and (|3.16p . If g satisfies (F2), by integrating (j3.3p with respect to time, we 
arrive at the conclusion. □ 

3.2 Higher-order energy estimates 

Denote 

A{t) = v\\Vv{t)f + \\M{t) - f{d{t))\?- (3.17) 
Lemma 3.2. The following inequality holds for smooth solutions {p, v, d) to problem (jl.ip - 

j^A{t) + Wp^m? + 7l|V(Ad(t) - f{d{t)))f 
< C{A\t) + A{t)) + C\\gf, Vt>0, (3.18) 

where C is a constant depending on \\vQ\\,\\dQ\\iii ,r], p, pjVjil. and also \\4>\\h^ (under (Fl)) 
l|g|lL2(o,+oo;L2) (under (F2)). 

Proof. Using equations ()l.ip - ()1.4p and the facts ()3.7p . ()3.8p . we compute that 

V d ..^ ..n fx 1 

— r = — / uAv-Vidx 

p\vt\'^dx— / p{v ■ Vv) ■ vtdx 
n Jn 
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-A / [{Ad- f{d)) ■Vd]-vtdx+ / pg-vtdx (3.19) 
Jn Jn 



and 



{Adt - f{d)dt) ■ {Ad - f{d))dx 
n 

' A{v ■ Vd) ■ {Ad - f{d))dx - -f\\V{Ad - 



- / f'{d)[-vVd + j{Ad-f{d))]{Ad-f{d))dx 
Jn 

- [{Av Vd){Ad - f{d))dx - f [{v ■ \7)Ad + 2VvV^d\ ■ {Ad - f{d))dx 
J Jn 

-7||V(Ad - f{d))f - 7 / f'{d){Ad - f{d)) ■ {Ad - f{d))dx. (3.20) 



n 



Adding p.l9p and ()3.20p together, we have 

~ (HlVt^f + \\Ad- f{d)f) + p\vt\''dx + 7||V(Ad - f{d))f 

= - [ p{v- Vv) ■ vtdx - [ [{Ad - f{d)) ■ Vd] ■ {Av + \vt)dx 
Jn Jn 

+ / pg-vtdx- / {vV)Ad- {Ad - f{d))dx 
Jn Jn 

-2 / {Vv\7'^d)-{Ad- f{d))dx--f [ f'{d){Ad - f{d)) ■ {Ad - f{d))dx 
n 
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^rn- (3.21) 



m=l 



Next, we estimate term by term. Using the fact (|3.7p . we rewrite (|1.2p as 

- lyAv + V ( P + ^ I Vdp + XF{d)] = -p{vt + vVv)- \{Ad - f{d))-Vd + pg. (3.22) 



2' ' '7 
By Lemma \2.1\ we get 

II^IIh^ < C{[[pvt[\ + [[p{vVM + [[{Ad-f{d))-Vd[\ + [[pg[\) 

< C{p-2 [[p-2vt[[ + p||^||l4 II Vi;||L4 + p||g||) + C||(Ad - f{d)) ■ Vd[[ 



< C{[[p-2vt[[ + C||7;|| 2 ||V7;||||^;||^, + ||g||) + C||(Ad - /(d)) • Vd|| 
1 

T 



< hMH^+C{[[p'^vt[\ + [\g[\ + [\{Ad- f{d))-Vd[\) + C[\Vv[\^ (3.23) 



namely, 

||7;||h2 < C{[[p'2vt[[ + ||g|| + ||(Ad - f{d)) ■ Vd[[) + C[[Vv[[\ (3.24) 
On the other hand, we have 

[[{Ad-f{d))-Vd[\ 

< ||Vd||L4||Ad-/(d)||L4 



< C{\\Adp\\Vdp + \\Vd\\)\\Ad- f{d)p\\V{Ad- fid))p 

< C{\\Ad-f{d)\\-2+l)\\Ad-f{d)\\-2\\V{Ad-f{d))f^. (3.25) 

Then we infer from p.24p . ()3.25p and the Young inequahty that 



^ hi A'' 



n 



'dx + 4 / p\v ■ Vupdx 

'dx + 4p\\v\\l4\\Vv\\l4 

'(ix + C7||z;||||Vw||(||Vf||||A?;|| + \\Vvf) 

'dx + C\\v\\ii2\\Vvf + C\\Vvf 

■dx + C\\p^vt\\\\Vvf + C(||g|| + \\Vvf)\\Vvf + C||V?;f 



2 



< 



+C{\\Ad - f{d)p + l)||Ad - /(d)||2 ||V(Ad - f{d))p\\Vv 
^ ^ pbtpdx + |i|V(Ad - f{d))f + C\\Vvf + C\\Vvf 
+C\\Ad - f{d)f + C\\Ad - f{d)t + CM\\ 



h < Ci\\Av\\ + p2\\p-2vt\\mAd-fid))-Vd\\ 

< C{\\p-2vt\\ + llgll + \\{Ad - fid)) ■ Vd\\ + \\Vvf)\\{Ad - fid)) ■ Vd\\ 

< Ci\\p^vt\\ + \\g\\ + \\Vvf) 

x(||Ad - fid)f^ + l)||Ad - /(d)||5||V(Ad - fid))\\h 
+Ci\\Ad - fid)f + \\Ad - fid)\\)\\ViAd - fid))\\ 

< y^p\v,\^dx + ^\\ViAd-fid))f 

+C\\Ad - /(d)r + C||Ad - fid)f + C\\Vv\\^ + C||gf. 

Concerning I^, it fohows that 

h<P^M\\p'^vt\\ < ^Wpvtf + 2p\\gf. 
Next, using ()3.14p . we get 

\\f'id)Vd\\ < CiMls + l)||Vd||L4 < C(l + \\Ad - /((i)ll^), 
which imphes that 

h < C||^;||L4||VAd||||A<i-/(d)||L4 

< C\\vf2\\Vv\\k\NiAd-fid))\\ + \\f'id)Vd\\)\\Ad-fid)\\-2\\ViAd-fid))\\ 

< C\\v\\-2\\Vvf^\\Ad-fid)f^\\ViAd-fid))\\l 
+C\\vf^\\Vv\\H\\^d - fid)f^ + \\Ad- /(d)||)||V(Ad - /(d))||l 

< |||V(Ad - fid))f + C\\Ad - fid)\f + C\\Ad - fid)f 

+C\\Vv\\^ + C\\\/v\\\ 
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From the elliptic estimate and ()3.14p . we have 



||d||H2 < C^IIAdll + lldll^a^^^j <C(||Ad-/(d)|| + ||/(d)|| + ||do||HO 

< Ci\\Ad-fid)\\ + l). (3.26) 
Then it follows from (j3.24p . (j3.25p . ()3.26p and Young inequality that 

h < C||V7;||L4||d||H2||Ad-/((i)||L4 

< C||V^||5||^|||,(||Ad-/((i)|| + l)||Ad-/(d)||^||V(Ad-/(d))||5 

< C\\Vv\\^^{\\p'2vtf2 + ||g||l + \\Vv\\)i\\Ad - fid)\\ + 1) 
x\\Ad-fid)\\^\\V{Ad-f{d))f2 
+C\\Vv\\-2md - f{d)\\i + l)||Ad - ||V(Ad - 

< P\vt\'dx + |||V(Ad - f{d))f + C||gf 
+C\\Vv\\^ + C\\Vvf + C\\Ad - + C\\Ad - f{d)\\\ 

Finally, for Iq, we have 

h < C||/'(a!)|lL4||A(i-/((i)||L4||Ad-/(d)|| 

< C{\\d\\ls + l)||V(Ad - /(d))||^!|Ad - 

< l\\V{Ad-fid))f + C\\Ad-f{d)f. 

Collecting the estimates for Ii, we infer from (|3.2ip that 

-A{t) < A f p\vt\-'dx-^\\V{Ad-f{d))f + C\\Vv\f + C\\Vv\f 



1 d 
Ydt 

+C\\Ad - f{d)f + C||Ad - + C||g||^ 



which yields our conclusion ()3.18p . □ 

Proposition 3.2. Under the assumptions of Theorem \l.l\ the following uniform estimates 
hold for any t > 0; 

\Ht)U. + \\d{t)\\^2<C, (3.27) 
sup / \\p^vt(r)f + \\V{Ad{T) - f{d{T)))fdT < C, (3.28) 



i>0 Jt 



where C is a constant depending on HwoHhii \\do\\u'^,r], p, p and also \\(P\\h'^ (under (Fl)) 
or ||g||L2(o,+oo;L2) (under (F2)). 

Proof. For both the cases (Fl) and (F2), (|3.15p implies that 

/ A{t)dt < +00. 
Jo 

The uniform bound ()3.27p follows from Lemma 13.11 the higher-order energy inequality 
()3.18p and the uniform Gronwall lemma [301 Lemma IILl.l]. Integrating ()3.18p from t to 
t + 1, we can conclude ([OH]) from ([3211) • □ 
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Denote the quantity 

B{t)= [ p\vt{t)\^dx+\\Vdt{t)f. (3.29) 
We can derive the following higher-order differential inequality: 

Lemma 3.3. // g satisfies (Fl), then following inequality holds for smooth solutions 
{p,v,d) to problem (|l.ip - ()1.6p ; 

^B{t) + i^\\Vvtf + j\\Adtf <C{B^{t) + B{t)) + C\\S/vf, Vt > 0, (3.30) 

where C is a constant depending on z^, ||uo||h1; MoIIh^j^jPjP o.i^d \\(j)\\jj2. On the other 
hand, if g satisfies (F2), then we have: 

j^B{t) + u\\Vvtf + j\\Adtf < C{B\t) + B{t)) + C(||g||2j, + ||g,f + ||Vt;||2), Vt > 0, 

(3.31) 

where C is a constant depending on i/, HvqUhIi MoIIh^j^/jPjP o^c? |lg||L2(Q__,_oc. . 

Proof. Taking temporal derivative of ()1.2p and taking inner product of the resultant 
with vt, then using the fact S/ ■ v = S/ ■ vt = and the transport equation (jl.ip . after 
integration by parts (keeping in mind that v\r = vt\r = 0), we get 



^ / pt\vt\'^dx-^ [ pv ■ V {\vt\'^)dx - [ p{vt ■ V)v ■ Vtdx 
Jn ^ Jn Jn 

- / pt{vt + v - Vv) ■ Vtdx - I VPfVtdx+ / {ptg + pgt) ■ vtdx 
Jn Jn Jn 

+2A / {Vdt Vd) : Wvtdx 
Jn 

- / pt\vt\^dx- / pt{v - Vv) ■ Vtdx - / p{vt - Vv) ■ Vtdx 
Jn Jn Jn 

+ / (Ptg + pgt) ■ Vtdx + 2A / (Vdt Vd) : Vvtdx 
Jn Jn 

- / pv ■ \/{\vt\'^)dx — / pv ■ V^v ■ Vv) ■ vt]dx — / p{vt ■ V v) ■ vtdx 
Jn Jn Jn 

+ pv- V(g • vt)dx + / pgf Vtdx + 2A / (Vdt Vd) : Vvtdx 
Jn Jn Jn 

6 

^ J^. (3.32) 



m=l 



First, we consider the case that g satisfies the assumption (Fl). In this case, J5 = 0. 
Using the uniform estimates in Proposition 13.21 we have 

Ji < Cp\\Vvt\\\\vth4v\\-i^4<C\\Vvt\\{\\Vvt\\^\\vt\\^ + \\vt\\) 
< ^W'^vtf + CWvtf. 

Then for J2, we infer from (fCTll . (l3:25]l and (13:2711 that 

J2 < P(l|Vf ||L4||f IlL^II^'tllL* + IbllLsll'^llHHktllL^ + ll''^llLs||Vt'||L4||Vft||) 
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< C[(||At>||||Vt;|| + ||V^;f )||V7;|| + IkHnHlVt^f 
+ \\Vvf{\\Avf2\\Vv\\^ + ||Vf|l)]||Vt;t|| 

< C\\Vv\\\\Vvt\\ + C\\v\\n2\\Vv\\\\Vvt\\ 

< C\\Vv\\\\Vvt\\ + C{\\p^vt\\ + ||g|| + ||Vt;f )||Vt;||||V7;t|| 
+C(||Aci-/(d)||5+l)||Ad-/(ci)||i||V(Ad-/(d))||i||V^;||||Vt;,|| 

< ^INvtf + C||V(Ad - fid))f\\Vvf + C\\phtf\\Vvf 

+ ||gf )||Vt;f . (3.33) 

When g = V(p, then the last term on the right-hand side of (|3.33|) is controlled by C|| Vu|p. 
Next, 



J4 = pv ■ V(V(/) • vt)dx 
Jn 

< P{\HL4vth4nH^ + l|V^t|||b|lL4||V0||L4^ 



< C\\Vv\\\\S/vt\\ 

Finally, for Jg, we have 

Je < C\\Vd\\L^\\Vdt\\\\Vvt\\ 

< C\\Vd\\H2\\Vd\\\\Vdt\\\\Vvt\\ 

< C{\\\7iAd-f{dm + l)\\Vdt\\\\Vvt\\ 

< ^llV^tf + CiWViAd - fid))f + l)\mf. 

Taking temporal derivative of (jl.4p and taking inner product of the resultant with 
—Adt, we obtain that 

i|||V*||^+.||Ad,|P 

= {vf Vd) • Adtdx + {v- Vdt) ■ Adtdx + 7 / f'{d)dt ■ Adtdx 
Jn Jn Jn 

:= J7 + J8 + J9, (3.34) 

where the right-hand side of p.34p can be estimated as follows 

J7 < \\vth4'^dh4Adt\\<C\\Vvt\\^\\vt\\'^\\Adt\\ 
< ^\\Vvtf + l\\Adt\\ + C\\vtf, 

Js < \\vh4Vdth4Adt\\<C{\\Vdtt^\\Adt\\^ + \\Vdt\\)\\Adt\\ 

7|| A J Il2 ^iiwj Il2 

u;||Loo||o.t||||^o.t|| ^ 

It follows from (fO]) that 



< ]:\\Adtf + C\\Vdtf, 
b 

J9 < 7\\f'{d)ho.\\dt\\\\Adt\\<'^\\Adtf + C\\Vdtf. 





Vdt = -V{v-Vd)+-f\/{Ad- f{d)). (3.35) 
13 



Using the uniform estimate (|3.27p and Sobolev embedding theorem, we have 

\\V{vVd)\\ < C\\Vv\\\\S/dUoo+C\\vh4V^d\\j^, 

< C||VHI||Vd|||,||Vd||Uc||V^||(||VAd||i||Ad||i + ||Ad||) 

< c||Vt;||||VAd||5 + (:;||Vf|| 

< C||Vt;||||V(Ad- +C7|IV7;|I, (3.36) 
which together with the equahty (j3.35p imphes that 

||V(Ad-/(d))|| < l\\Vdt\\ + -\\V{vVd)\\ 

7 7 

< illVdtll + C||Vt;||||V(Ad - f{d))f2 + C\\Vv\\ 
7 

< -||Vdi|| + i||V(Ad-/(d))||+C||Vt;||. (3.37) 
7 2 

As a result, we get 

llVdfll < ||V(^;- Vd) 11+711 V(A(i-/((i)) II 

< C||Vz;||||V(Ad - /(d))||i + Clival + 7l|V(Ad - /(d))|| 

< (l + 7)||V(Ad-/(d))||+C||V7;||. (3.38) 

Collecting the above estimates for Ji,...,Jq and using the relations ()3.37p . ()3.38p . we 
conclude from (j3.32p and (j3.34p that 

±B{t) + u\\Vvtf + j\mf 

< C\\Vvf{\\phtf + \\Vdtf) + C\\Vdtf + C{\\vtf + llVdtf ) + C7||V^;f , 

which yields our conclusion ()3.30p . 

Now if g satisfies (F2), we only have to re-estimate the terms J2, J4 and J5 using the 
uniform estimate (|3.27p . It follows from (|3.33p that 

J2 < ^WVvtf + C\\V{Ad - f{d))f\\Vvf + C\\pht\\'\\Vvf + C{1 + ||gf )||Vt;f 

< ^llV^tf + C||V(Ad - f{d))f\\Vvf + C\\p-2vtf\\Vvf + C\\gf + C\\Vvf. 

Moreover, 

J4 < ^(lbl|L4|kt|lL4||g||Hi + l|Vft||||t;||L4||g||L4) 

< C^llg|lHi||||V^|l||Vt;i|| 

< ^llV^tf + C^llg|l^i, 

J5 < fclllKII<^l|Vt;if + C||g,f. 

By the above estimates, we easily arrive at ()3.3ip . □ 

Proposition 3.3. Under the assumptions of Theorem \l.l\ the following uniform estimates 
hold for any t > 0, 

\\vt{t)\\ + \\dt{t)\\^.<C, (3.39) 
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vm^2 + \\dm^s<c, (3.40) 
t+i 



sup 

t>0 Jt 



sup 

t>0 Jt 



/t+i 

i\\VvtiT)f + \\AdtiT)f)dT<C, (3.41) 

'^*^'lkWllw2>.^-^<'^' VgG(l,+oo). (3.42) 
Moreover, for any T > 0, 

\\p{t)\\w^,, < e^^WpoWwi,,, Vt G [0,r]. (3.43) 

The constant C in the above estimates depends on IjfollH^) WdoW'R^^'H^ P^P ^^^^ W'PWh'^ 
(under (Fl)) or ||g||L2(o .(.j^. jji), ||gt||L2(o,+oo;L2) (under (F2)), hut it is independent of t. 

Proof. It follows from the definition of B{t) and ([338]) that 

rt+l 



/t+i 
B{T)dT < C, (3.44) 



sup 

t>0 Jt 

where C is a constant depending on z/, ||t;o||h1) Mo||h2)^)P)P and also ||(/>||//i (under (Fl)) 
l|g|lL2(o,+oo:L2) (under (F2)). Applying the uniform Gronwall lemma again, we infer 
from (fSrHD]) or (fOTD . and (imH that 

B{t) <C, Vt > 0, (3.45) 

where C is a constant depending on u, H^'oIIh^) \\do\\u^,r], p, p, and also |[</>||/i-2 (under (Fl)) 
l|g|lL2(o,+oo;Hi)) l|gt IIl2(o,+oo;L2) (under (F2)). This yields the estimate (|3.39p . Inte- 
grating (|3.30p or (|3.3ip from t to t + 1, we conclude (|3.4ip . RecaUing ()3.22p . we deduce 
from Lemma 12.11 and the estimates ()3.27p , ()3.37p that 

\\v\\n^ < C{\\pvt\\ + \\p{vVv)\\ + \\{Ad- f{d))-Vd\\ + \\pg\\) 

< Cp{\\vt\\ + |b||L4||Vt;||L4 + llgll) + C||Vd||L4||Ad- /(d)||L4 

< C{\\vt\\ + \\Vv\\\\v\\l,\\Vvf2 + \\g\\) 
+C||d||H2(||V(Ad-/(d))||i||Ad-/(d)||l) 

< \\Hm+C{\\vt\\ + \\ydt\\ + \\Vv\\ + \\g\\). (3.46) 

Using ()3.39p and the facts that either under (Fl) ||g|| < or under (F2) ||g|| is 

bounded by a constant depending on ||g||L2(o,+oo;L2) and ||gt||L2(o,+oo;L2)> we get the uni- 
form estimate for ||f ||h2- 

Next, by the elliptic estimate, we infer from (|3.35p that 

||Vd||H2 < C (llVdill + \\Viv . Vd)|| + \\nd)Vd\\ + l|Vdo||^3^^^ + ||Vd|l) 

< C + C||V7;||L4||Vd||L4 +C||t;||L-||d||H2 

< C7 + C||^||H2||d||H2. (3.47) 

Combining ()3.46p and ()3.47p . we obtain the uniform estimate ()3.40p . In a similar manner 
to ^M\\ . for q G (1,-Foo), using (fOO]) . we get 

\\v\\w^,, < C{\\pvt\\-L, + Mv ■ Vv)\\-L, + \\{Ad - f{d)) ■ Vc^IIl, + \\pghi) 
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< CpiWvtU, + \\v\\^. \\Vv\\^. + ||g||L.) + C||Vd||^| ||Ad - f{d)\\^. 

< C{\\Vvt\\ + \\v\\^2 + ||V(Ad - f{d))\\ + llgllL,) 

< C{\\Vvt\\ + \\V{I^d- f{d))\\ + M\^. + \\Vv\\). (3.48) 

Thus, (|3.42p is a consequence of (|3.48|) . (|3.28p and (j3.4ip . where the bound either depends 
on \\<i)\\m (under (Fl)) or ||g||L2(o,+oo;Hi) (under (F2)). Finally, for any T > 0, 

^ ||Vf(T)|lLoodr < C \\v{T)\\^2,rdT<CT-2 (^j^ HT)f^2.rdTy 

< CT, Vt G [0,r], for some r > 2. 
Then we infer from (|2.20p that (|3.43p holds. The proof is complete. □ 

3.3 Proof of Theorem [TTT] 

Under the assumptions in Theorem 11.11 we can derive estimates for the approximate 
solutions to problem ()l.ip - (|1.6p {pm,Vm,dm) as in Propositions 13. 1| 13.21 and 13.31 which 
are independent of the parameter m. By extracting a subsequence and passing to limit 
as m — >■ +00, we can obtain a global strong solution {p,v,d) to problem ()l.ip - ()1.6p in a 
standard way such that 

pGL'^{[0,T],W'''{n)), 

V eL'^i[0,T];U\n)nV)nL\0,T;W^^^), t>t G L2(0,r;V), gG(l,+oo), 
d G L°°([0,r];H3(f))), dt G L°^{[0,T]-Ul{n)) n L\0,T;U^), 
< p< p{x,t) < p, \d{x,t)\<l, V(x,t) G X [0,T]. 

Due to Lemma 12.21 we infer that p G C{[0,T],W^'^ {Q)). It is not difficult to see from 
(f3:i0]l and (IHIiTTl that d G L'^{0,T;W^), which together with dt G L'^{0,T;ll^) yields that 
d G C([0, T]; H'^). By the regularity of the Stokes operator, 

II^IIh3 < C{\\pvt\\ni + \\p{v ■ Vv)\\ni + \\{Ad - f{d)) ■ Vd\\^, + \\pg\\ni) 
< C'/odlutllni + \\v ■ VvIIhi + ||g||Hi) 

+C\\Vp\\-^,,{\\vt\\-L, + 11^; • VvW-L, + ||g||LO + C||V(i||L-||Ad - f{d)\\ui 

+C\\V'^d\\j^g>\\Ad - f{d)\\i^q, for some 2 < < +00, - + ^ = i 

q q 2 

and estimates (lOOll . (lOTT) . (103]) . we infer that v G L'^{0,T;1l'^). This fact and vt G 
L2(0,r; V) yields the continuity v G C([0,r],H2). 

Finally, we briefly show that the strong solution is indeed unique. Suppose {p, v, d) 
and {p,v,d) are two strong solutions corresponding to the same initial data {po,vo,dQ). 
One can easily check that all the computations in [3,, Section 4] can be verified due to the 
regularity of the two strong solutions. Denote 6p = p — p, 6v = v — v and dd = d — d. We 
have for any t G [0,r] (cf. e.g., [31 (4.46)]): 

l\Mt)f + p{t)\6v{t)\^dx + \\\V{6d{t))f 

< - [ \\Vi5v)fdt- [ \\A{dd)fdt+ [ I {5p){5v)Vpdxdt 
Jo Jo Jo Jn 
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pVv\5v\^dxdt+ [ [ 
Jo Jn 



{5p){5v){vt + V ■ Vv)dxdt 



'0 Jn Jo Jn 

+ [ [ v{5d)A{5d)dxdt - f f i6v)V{5d)Addxdt 
Jo Jn Jo Jn 



+ 



/ A{5d)ifid) - f{d))dxdt. 
Jn 



Using the Holder inequality, Young inequality and the estimates ()3.39p , ()3.40p , ()3.43p , we 
easily get 

¥p{t)f+ f pmv{t)\^dx+\\V{5d{t))f <Ct [\\\5pf + \\Svf + \\Vi6d)f)dt. 
Jn Jo 

Due to the positivity of the density p > p > 0, the uniqueness follows from the above 
estimate and the Gronwall inequality. The proof is complete. □ 

4 Long-time behavior 

In this section, we study the long-time behavior of global strong solutions to problem 

dLiD-dini). 

Proposition 4.1. Under the assumptions in Theorem \l.ll either g satisfies (Fl) or (F2), 
the global strong solution to (jl.ip - ()1.6p has the following decay property: 

hm (|l^;(t)||Hi + ||Ad(t)-/(d(t))||)=0. (4.1) 

t— J- + 00 

Proof. For both cases (Fl) and (F2), (fXT^ implies that A{t)dt < +oo. Then it 
follows from (|3.18p and [35, Lemma 6.2.1] that limt^+oo^(i) = 0, which together with 
the Poincare inequality yields ()4.ip . □ 

Proposition 4.2. Under the assumptions in Theorem \l.l[ either g satisfies (Fl) or (F2), 
the global strong solution to (jl.ip - (|1.6p has the following decay property: 

hm {\\vt{t)\\^, + \\V{Ad{t) - f{d{t)))\\ + lldtllHi) = 0. (4.2) 

Proof. First, we look at the case that g satisfies (F2). We have already proved that A{t) 
is bounded for all time and A{t) G L^(0,-|-oo). Therefore, integrating (|3.18p from to 
+00 with respect to time, we infer from (F2) that 

+00 

{\\p--vm\^ +i\\^{^d(t) - f{dm\?)dt 







< ^(0) + C f sup A{t) + l] r A{t)dt + C [ 
\t>o J Jo Jo 

Recalling (j3.38p . we obtain 



-00 



\g{t)fdt<+oo. (4.3) 







B{t)dt < +00. (4.4) 
Then using (|3.3ip . assumption (F2) and [351 Lemma 6.2.1], we conclude that 

lim B{t) = 0, (4.5) 

t— 7- + 00 
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which combined with Poincare inequaUty and (|3.37p . (|3.38p yields (|4.2p . 

Next, we deal with the case when g satisfies (Fl). We only need to show that ()4.4p 
still holds in this case. For this purpose, we re-estimate the term I2, I3, I4 and Is in (|3.2ip 
using the higher-order estimate (j3.40p instead of the lower-order one (j3.27p . As in |34J , we 
infer from the transport equation (|l.ip and integration by parts that 

Is = — /" pVcj) ■ vdx — [ ptVcj) ■ vdx 
dt Jn Jn 

= — / pV(f> ■ vdx — pv ■ V(V(/) • v)dx 
dt Jn Jn 
d f 

d f 

< — / pV(l) ■ vdx + C{\\v\\\\Vv\\ + \\v\\^\\Vi 
dt Jn 

< ^ [ pV(l)-vdx + C\\Vvf. 
dt Jn 

Next, since Ad — f{d)\r = 0, we obtain that 



'v\\ 2 



h + h + h = -A 



f [{Ad - f{d)) ■ \7d] ■ vtdx - [ A{v Vd) ■ {Ad - f{d))dx 
Jn Jn 

X [ [{Ad - f{d)) ■ Vd] ■ Vtdx + [ \/{v Vd) ■ \/{Ad - f{d))dx 
Jn Jn 



= -X [{Ad - f{d)) ■ Vd] ■ Vtdx + / V kViV idjV k{Adj - fj{d))dx 
Jn Jn 

+ / ViVkVidjVk{Adj - fj{d))dx 
Jn 

:= + + (4.6) 

where 

I2 < Cp^\\p-2vt\\\\{Ad- f{d))-Vd\\ 

< C\\p'^vt\\{\\Ad - + l)||Ad - ||V(Ad - 

< ^ p\vt\^dx + |||V(Ad - f{d))f + C\\Ad - + C||Ad - f{d)f, 

+ < \\Vv\\\\Vdhoo\\V{Ad-f{dm + \\vh4V^dh4V{Ad-fidm 

< C||V7;||||d||H3||V(Ad-/(<i))|| 

< l\\V{Ad-f{d))f + C\\Vvf, 

Replacing the original estimates for l2,---i H by -^2v) -^5; we arrive at the following in- 
equality 



^iA{t)-2 / pV<l)-vdx]+\\p^vt{t)f + ^\\V{Ad{t)-f{d{t)))f 



dt 

< C{A'{t) + A{t)), (4.7) 

where C is a constant depending on \\vo\\, \\do\\iii,r], p, p,i^,Q and also ||0||j:/2. Integrating 
()4.7p with respect to time from to +00, we deduce that 

\p'2vt{t)f +j\\V{Ad{t) - f{d{t)))fdt 



18 



< ^(0) - 2 / poV0 • vodx + 2 pV(t)vdx 
Jn Jn 

+C (sup A{t) + l] r A{t)dt 
\t>o J Jo 



< +00, 

which again imphes ()4.4p . Using the same argument as for the previous case, we conclude 
the decay property ()4.2p . The proof is complete. □ 



4.1 Proof of Theorem Q 

According to Propositions 14.11 and 14.21 it remains to show the convergence for the density 
function p and the director vector d. 

As a direct consequence of the uniform-in-time estimate (j3.12p and weak compactness 
of bounded sets in L'^ (1 < q < +oo), we know that for any sequence {ti} +cx3, there is 
a subsequence {t^} +oo such that p(t'^ weakly converge to a certain p^o in L'^ . On the 
other hand, due to Lemma YTHX \p(t'^)\Lq = ||poo||l9 = ||/Oo||l9; since L"? (1 < g < +oo) is 
a uniformly convex Banach space, we can conclude that p{t'j) actually strongly converge 
to /9oo in L'^ (cf. |351 Lemma 3.1.6]), namely, (|1.14p holds true. 

The uniform- in-time estimate ()3.40p yields that for any sequence {ti} +oo, there is 
a subsequence {t^} +cx3 such that d(t'^ strongly converge to a certain d^o in H^. Then 
by ()4.ip . we see that 

< II - Adoo + /(doo)|| 

< II - Ad(t^ + fWM + - Adooll + ll/KtD) - /((ioo)ll 

< II - Ad(t^ + /(d(t^))|| + C\W,) - doollH^ 

0, as t ■ /- +00. 

Thus, doo satisfies the stationary problem ()1.16p . The convergence follows from ()4.2p 
and the fact that 

||VA(d(tO - d^)\\ < \\V{-Ad{t',) + fidmn + l|V(/KtD) - /(doo))|| 
< ||V(-A(i(tO + fid{t'M\ + C\\d{t'^ - doollH^ 
0, as +00. 

The proof of Theorem 11.21 is complete. 

Remark 4.1. (1) When the external force g is a gradient field independent of time, our 
results show that the velocity field and its time derivative will converge to zero as time 
goes to infinity, this coincides with the result for the 2D density- dependent incompressible 
Navier-Stokes equations J3^. 

(2) For the density function p, we are only able to obtain some partial results that it will 
sequentially converge to a certain function in L'^ norm (differential sequences may have 
different limit points). One possible sufficient condition for the convergence in time is that 
^ L^(0,+oo). On the other hand, we observe that for any function ip € Hq{^), 

pt^pdx = — I V ■ {pv)ipdx = / pv ■ Vipdx < p\\v\\\\V(p\\, 
n Jn Jn 
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which yields WptWn-'^ ^ C'Ibll- Hence, if we can prove \\v\\ € L^(0, +00), then we get the 
convergence of p as time tends to infinity in the space and thus in L'^ due to the 

uniqueness of limit. However, we do not know this -integrability condition on v from 
the above proof. 

(3) For the director vector d, we are able to show the decay of its time derivative and 
sequential convergence of itself to a steady state that is a solution to the stationary problem 
()1.16p . If we know that problem ()1.16p admits a unique solution, then d will converge to it 
as time goes to infinity. However, in general we cannot expect the uniqueness of solutions 
to the stationary Ginzburg-Landau equation ()1.16p . unless, for instance, the parameter r] 
in f{d) is sufficiently large. 

4.2 Proof of Theorem [TTSl 

In this subsection, we provide the proof for Theorem II .31 Now the external force g is an 
asymptotically autonomous one such that it satisfies (F2) and (|1.17p . Different from the 
previous case with time-independent force, we shall see that one can obtain convergence 
for the density function and director vector. The proof is based on an appropriate gener- 
alization of the Lojasiewicz-Simon approach for gradient-like systems (cf. e.g., [T2|[T3]). 
Denote 

E{d) = \\\Vdf + jj{d)dx. (4.8) 

It is straightforward to check that any solution to the stationary problem (|1.16p is a critical 
point of the energy functional E{d), and conversely, the critical point of E{d) is a solution 
to ()1.16p (cf. [33J). Besides, regularity of solutions to (|1.16p has been shown in [20] such 
that d is smooth on Q provided that do is smooth on F. Below we shall make use the 
following Lojasiewicz-Simon type inequality (cf. [T31I33]): 

Lemma 4.1. Let ifj be a critical point of E[d). There exist constants € (0, ^) and /? > 
depending on such that for any d € H^(J7) satisfying d\r = do{x) and \\d — V'IIh^ < P> 
there holds 

\\-Ad + f{d)\\ > \E{d) - E{i;)\^-'^ . (4.9) 

The uniform estimate (|3.40p yields that the set of all limit points of the trajectory 
{d(t) : t > 0} denoted by uj{d) := ns>o {^(0 ^ : t > s} is non-empty and compact 

in H^(r2). Besides, similar to Section 4.1, we can see that each element in io{d) is a solution 
to problem ()1.16p and thus is a critical point of E{d). We infer from p.3p that 

rti 

£{ti) -£{t2) <C Wgfdt, Vti>t2>0. (4.10) 

Jt2 

Since ||g|| G 1/^(0, +cc), it follows that S{t) converges to a certain constant £qo as t — )■ -|-oo. 
Recalling the convergence of velocity field ()4.ip . we deduce that E{d) is indeed a constant 
on u{d) such that E{ip) = X~^£oo, for all G uj{d). Due to Lemma ini for every G u:{d), 
there exist some and 9^ G (0, |) that may depend on tp such that the inequality (|4.9p 
holds for d G B/s^iiP) := {d G : d\r = do, \\d - ^pWu^ < /3^} and \E{d) - E{'tP)\ < 1. 

The union of balls {B|3^{^p) : ip G oj{d)} forms an open cover of ijj{d) and due to the 
compactness of uj{d), we can find a finite sub-cover {B/3.(-0j) : i = 1,2, ...,m} where the 
constants /3, 9 corresponding to ipi in Lemma l4.ll are indexed by i. From the definition of 
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ijj{d), we know that there exist a sufficient large to such that d{t) £ U := U™ ^^S/?, (V'j) for 
t > to. Taking 6 = min^jej € (0, i), we get for ah t > to 

II - Ad{t) + fidmi > \E{d{t)) - X-'£oo\'-'. (4.11) 

Denote 

y{tf = -\\\7v{t)f + Xj\\Ad{t) - f{d{t))f, m= / Mr)\\'dT. 

J t 

Assumption (|1.17|) imphes that z{t) < C(l + 1)~^-^~^^^ for t > 0. Then by the basic energy 
inequahty ()3.3p . we get 

£{t)-£oo> y{rfdT--^z{t)> y[rfdT-C{l + t)'^^+^\ (4.12) 

On the other hand, using (j4.1ip . the uniform estimates (|3.27p and the fact < 2, we 
obtain that 

|^:(t)-^oo| < \\\vf + \\E{d)-X-^£oo\ 

< C\\Vv\\^ +\\\- Ad{t) + f{d{t))\\^ 

< cy{t)^, yt>to. (4.13) 



C = min <; e, — - 1 G f 0, - 

^ ^ ' 2(1 + J V 2 



Take 

It is easy to check that 

/•+00 r+oo 

/ (i + ^)-2(i+?)(i-C)rf^< / (l + r)-(2+«)dr < (l + t)-(i+«), Vt>0. (4.14) 
Jt Jt 

We now set 

Z(t) = 3^(t) + (l + t)-(i+«)(i-^). 
Since hmt_,+oo 3^(t) = (see (gT])), it fohows from KWi . ()i33D and ()mD that 

/ + 00 ^ J 

zirfdT < cy{t)T=e +c{i + t)-^^+^^ <cy{t)— +c{i + ty'^^+^^ 

< CZ{t)^, yt>to. (4.15) 
Recall the following result (cf. [12^ Lemma 4.1] or ^ Lemma 7.1]) 

Lemma 4.2. Let ( E (0, ^). Assume that Z > be a measurable function on (0, +oo), 
Z £ L^(]R+) and there exist C > and to >0 such that 



P + OO ^ 

J Z{TfdT<CZ{t)~, for a.e. t > to- 



Then Z € L^{to, +oo). 
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We conclude from ()4.15p and Lemma 14.21 that 

r+oo 

/ Z{t) < +00. 

Jtn 



Since > 0, it holds 

»+oo 



r+oo r+oo 

/ (1 + < / (1 + t)-2(2+Odt < 2(1 + to)-^^+^^ < +00, for to > 0, 

Jtn Jtn 



'to -Jtn 
which implies that 



to 



\Vv{t)\\ + \\Adit) - fidit))\\dt < +00. (4.16) 

On the other hand, it follows from equation (|1.3p that 

\\dt\\ < C(i|^;||L4||V(i||L4 + II - Ad + f{d)\\) < C{\\Vv\\ + || - Ad + f{d)\\). (4.17) 
As a consequence, 

;>oo 

/ \\dt{t)\\dt < +00, (4.18) 
J to 

which easily implies that as t — > +00, d{t) converges strongly in By compactness 

of d{t) in H2(r2), we deduce that 

lim \\d{t) - dooW-n^ =0, (4.19) 

i— i>+oo 

where doo is a solution to problem (I1.16p . Recalling the uniform estimate (j3.40p . it holds 
llVAd- VAdooll < \\V{Ad-Ad^-f{d) + f{d^))\\ + \\V{fid)-f{doo))\\ 

< |lV(Ad-/(d))|| +C||d-(ioo||H2- 

The above estimate together with ()4.2p and ()4.19p yields 

lim \\d{t) - doollna = 0. (4.20) 

t-++oo 

It follows from Remark l4.ll and ()4.16p that ||/5t||/^-i G L^(to,+oo). Thus p{t) converges 
strongly in as t — > +00. By an argument similar to that in the proof of Theorem II. 2 ^ 
we conclude that 

lim ||p(t) -PooIIl? =0, gG(l,+oo). (4.21) 

t— >-+oo 

Next, we prove the convergence rate. Denote 

i2 ^2 r+°° 



fil ^2 f+OO 

}C{t) = S{t)-S^ + ^ l|g(r)fdr. 



It follows from the basic energy inequality (|3.3p that 

^K{t)+y{tf <Q. (4.22) 

Thus, K,{t) is decreasing on [0, +00) and K.{t) ^ as t — )■ +00. Recalling the definition of 
to, for t > to, we deduce from ()1.17p and ()4.13p that 

^(t)2(i-0 < C3^(t)^+C7(l + t)-2(i-^)(i+«) 
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< C3^(t)2 + C(l + i)-2(l-C){l+0 

< -C7^/C(t) + C(l + t)-2(i-?)(i+«), (4.23) 

where we have used the fact that '^^-^iP > 2. It foUows from the ordinary differential 
inequahty (|4.23|) and [21 Lemma 2.6] that 

IC{t) <C{l + ty\ Vt>to, (4.24) 

with the exponent given by 

L = min I , 1 + O = niin | , 1 + O ■ 



1-2C J tl-26' 

We infer from (|4.22p that for any t > Iq 



2t , / r'it \\ 



y{T)dT<t2( y{TydT] <ct^ic{t)2 <c{i + t) 



2 



where 



L-1 { e i\ 

mm < -, — > > 0. 

1 - 26*' 2 



It holds 



/ y{T)dT<y2 y{T)dT <cy2(.'^^t)-^ <c{i+ty, yt>to. (4.25) 

Then by (jHZD, we get 

I'+OD r+oo 

/ \\dt{T)\\dT< / yiT)dT <C{l+t)-^, Vt>to- 

Jt Jt 

which together with (j3.27p yields the convergence rate of d in 

\\d{t) -d^\\<C{l + t)-'', Vt>0. (4.26) 
Besides, using the Poincare inequality, we infer from Remark 14. II and (|4.25p that 

r+oo r+oo r+oo 

/ \\pt{r)\\H-^ < C \\vir)\\dT<C \\y{T)\\dr <Cil + t)-^, > 0, 

Jt Jt Jt 

which implies 

llp(i)-Pook-i <C(l + ^)-^ vt>o. 

Taking advantage of the lower-order convergence rate of the director d ()4.26p , we are 
able to obtain decay estimate for v as well as higher-order convergence rate on d. For this 
purpose, we make use the fact — Adoo + /(^oo) = and test ()1.4p by X{—A{d — doo) + 
f{d) - /(doo) + d- doo) to get 

■\\V{d - doo)f + ttM - dooW^ + A / F{d) - F{doo) - /(doo) • {d - doo)dx 



dt\2 



+A7II - Ad+/(d)f + A7||V((i-doo)f + A f {vS7d)-Addx 

Jn 
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= -X (v Vd) ■ (d - doo)dx - A7 / (fid) - fidoo)) ■ (d - doo)dx. (4.27) 
Jn Jo. 

Using the uniform estimate (|3.40p and Sobolev embeddings, we can estimate the right- 
hand side of (|4.27p as fohows 

R.H.S of (I07D < ^\\Vv\\^ + C\\d - doof . (4.28) 
Adding dOT]) with (I23D and ([23]), we infer from ([i:^ that 

j^Q{t) + CM{t) + l|V(d - doo)f ) < C2{\\d - doof + ||gf ), (4.29) 



where 



+A / F{d{t)) - F{doc) - fidoo) ■ {d{t) - doo)dx. 
Jn 

By the uniform estimate ()3.40p and the Taylor's formula, it is easy to see that 



= ll^P(tMt)\'dx + ^\\V{d{t)-doo)f + ^\m-doof 



Fid) - Fidoo) - fidoo)id - doo)dx 



n 



<Cs\\d-d 



2 

00 II ) 



which implies 



Qit) + C^Mt) -doo\\^>l [ pit)\vit)\^dx + ^i|V(d(t) - doo)\\^. (4.30) 

Since A(t) is uniformly bounded in time, ()3.18p can be rewritten as 

jAit)<C^Ait)^CM?. (4.31) 
We deduce from (IMD . M and (lOB that 

^Qi(t) + C6Qi(t) < Crdld - doof + ||g||'), (4.32) 
where Q\(t) = Qit) + j0^Ait). As a consequence, 

Qiit) < e-^«*(^Qi(0) + Cr j|*e^«^(||d(r)-doof + ||g(r)f)(iT) 

< Qi(0)e-^«* + Ce-%^ [\\\d{r)-doo\\^ + \\giT)f)dT 

Jo 

+c( sup ||d(s) -dooin e-^«* f e^^^dT + C f^^ \\git)fdT 
\se[{A J Jl Jl 



Cf:t 



< Qi(0)e-^«* + e-'^ ( C J\l + T)-^''dT + 

< (7(1 + *)"^", Vt>0. (4.33) 
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Recalling the definitions of A{t), Qi{t), Q{t), we conclude from (On|) . K26h and 

that 



\Ht)\\m + - ^^ooIIhi + \\^d{t) - f{dm\ < C(l + 0"", Vt > 0. (4.34) 
By the elliptic estimate 

||d-doo||H2 < C\\A{d-d^)\\<C\\Ad-f{d)\\+C\\f{d)-f{doo)\\ 

< C\\Ad- fid)\\ + \\d-doo\\m, (4.35) 

we get 

\\d{t) -d^\\ii2 KCil + ty, yt>0. (4.36) 
The proof of Theorem 11.31 is complete. 

Remark 4.2. If the velocity field v decays fast enough, we can obtain uniform-in-time 
W^'"^ -estimate for p {1 < r < +oo). For any q £ (1, +oo), we infer from (|3.48p . (j4.4p and 
(F2) that ||t;(t) ||-vv^2,9 € L^(0,+cx3). //k > |, or in other words, ^ > 3 and 6 S (|, ^), we 
just take q € (Ifrfi+oo) such that ^^3^!^^'" > 1- Then it follows from (|4.34|) that 

hoo 

\\Vv{t)\\l,oodt 
f + 9 „„2 

< c^/ wvmcxi'W'^vm^dt 



3i;-4 , 
, , 2,-4 \ 4(F^ / r+°° ^ \ 4(9-1) 

< C7f / \\Vv(t)\\—^dt] i / |b(t)||^2,,di 

3q-4 



< Ci {l + t)~^^dt\ 

< c. 

Therefore, by (j2.20p . we obtain that \\p{t)\\w^''- ^ C'UpollvFi-'' foi"t > 0. 
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